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ABSTRACT 

A multi-epoch, multi-instrument analysis of the Seyfert 1 galaxy HE 0436^717 is conducted 
using optical to X-ray data from XMM-Newton and Swift (including the BAT). Fitting of the 
UV-to-X-ray spectral energy distribution shows little evidence of extinction and the X-ray 
spectral analysis does not confirm previous reports of deep absorption edges from O Vlll. 
HE 0436^717 is a “bare” Seyfert with negligible line-of-sight absorption making it ideal to 
study the central X-ray emitting region. Three scenarios were considered to describe the X-ray 
data: partial covering absorption, blurred reflection, and soft Comptonization. All three inter¬ 
pretations describe the 0.5-10.0 keV spectra well. Extrapolating the models to 100 keV results 
in poorer fits for the the partial covering model. When also considering the rapid variability 
during one of the XMM-Newton observations, the blurred reflection model appears to describe 
all the observations in the most self-consistent manner. If adopted, the blurred reflection model 
requires a very low iron abundance in HE 0436^717. We consider the possibilities that this 
is an artifact of the fitting process, but it appears possible that it is intrinsic to the object. 

Key words: X-ray: galaxies - galaxies: active - galaxies: nuclei - galaxies: Seyfert - 
galaxies: individual: HE 0436-4717 


1 INTRODUCTION 

Seyfert 1 galaxies are a subclass of active galactic nuclei (AGN) 
that provide a preferred view of the central engine. Seyfert 1 galax¬ 
ies that appear to have little to no intrinsic absorption, so-called 
“bare” Seyfert Is, provide an especially unobscured line of sight to 
the innermost regions close to the supermassive black hole itself. 
Such AGN have been noted before (e.g. Vaughan et al. 2004, Em- 
manoulopoulos et al. 2011) and recently have become the subject 
of larger sample studies (Patrick et al. 2011, Jin et al. 2012, Walton 
etal. 2013). 

HE 0436-4717 appears at first to be a rather typical Seyfert 
1 galaxy at z = 0.053 first observed by the Einstein Slew Survey 
in 1992 (ESS, Elvis et al. 1992). It has been included in numer¬ 
ous surveys and catalogues since the ESS including: the ROSAT 
All-Sky Survey (Grupe et al. 1998), the Hamburg/ESO Survey for 
bright QSOs (Wisotzki et al. 2000), and the Palermo Swift-BAT 
Hard X-ray Catalogue (Cusumano et al. 2010). Despite its frequent 
inclusion in such wide-field studies, HE 0436-4717 has been the 
subject of deep X-ray analysis only once (Wang et al. 1998). 

Halpem & Marshall (1996) monitored the object in the ex¬ 
treme UV with the EUVE satellite for 20 continuous days and were 
able to produce a long light curve. The timing analysis of the light 
curve revealed possible flaring with a period of about 0.9 days. 
In addition, the authors also analyzed spectra from EUVE's short- 
wavelength spectrometer and found no emission line features in the 
70 - 110 A range. The column density listed in that work is compa¬ 


rable to Galactic extinction alone, indicating no intrinsic extinction 
in the object was observed. 

A detailed X-ray analysis of this object was done by Wang 
et al. (1998) using data from A5CA and ROSAT observations. The 
authors found the 0.4 - 10.0 keV spectrum was well modelled by 
a single power law with slope of T = 2.65±0.20 or by a sum of a 
r = 2.15±0.04 power law and a weak black body of temperature 
29.0±2.0 eV. The authors describe an Fe Kq emission feature with 
an equivalent width of around 430 eV and two absorption edges: 
one at £ ~ 0.28 keV with an optical depth of r ~ 3.46 measured 
from the ROSAT spectrum, and one at £ ~ 0.83 keV with r ~ 0.26 
from the ASCA spectra - the latter of which was attributed to a 
combination of O VII and O VIII edges. A comparison between 
ASCA observations 4 months apart saw a 47% flux change, however 
all other spectral parameters remained constant. 

In this study, an analysis of HE 0436-4717 multi-epoch data 
from XMM-Newton and Swift observatories spanning from 2002 
- 2013 is presented. Chandra also serendipitously observed this 
object between 2000 - 2006, however these data were determined 
to be piled up and were not used in the analysis (see Section |2]2](. 
The study of XMM-Newton and Swift data presented here is the 
deepest investigation of this object to date using the most current 
instruments and analysis tools. 

Following this introduction. Section describes the data pro¬ 
cessing for each observation. An initial X-ray spectral analysis of 
the XMM-Newton EPIC data is described in Sectionj^ including a 
search for absorption features seen previously in this object. The 
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Table 1. HE 0436-4717 Data Log. 


(1) 

Telescope 

(2) 

ObslD 

(3) 

Instrument 

(4) 

Filter 

(5) 

Start Date 

(6) 

Duration (s) 

(7) 

GTI (s) 

XMM-Newton 

0112320201 

pn 

thin 

09/10/2002 

68177 

67210 



MOS 1 

thin 


69167 

68270 



MOS 2 

thin 


69167 

68380 



OM 

UVWl 


34996* 

N/A 




V 


25000* 

N/A 


0603460101 

MOS 1 

thin 

12/15/2009 

129141 

126900 



MOS 2 

thin 


129146 

123600 

Swift 

00035763001 

XRT 

none 

12/07/2007 

10090 

9732 



UVOT 

V 


838 

838 




B 


840 

840 




U 


840 

840 




UVWl 


1681 

1681 




UVM2 


2286 

2286 




UVW2 


3369 

3369 


70-mo survey 

BAT 

none 

09/21/2010 

N/A 

N/A 


*lndividual filters exposed for a range of 1998 - 5001 s each with a sum of 10 exposures for UVWl and 5 exposures for V. 


spectral analysis was extended into the UV band in Section |4l 
which includes spectral energy distribution (SED) modelling using 
Swift data (UVOT, XRT, and BAT). SectionOdescribes an in-depth 
X-ray analysis of the XMM-Newton MOS spectra considering ab¬ 
sorption (Sec. HU, reflection (Sec. [Q l. and soft Comptonization 
(Sec. I5.4t scenarios. Section details the timing analysis of the 
XMM-Newton light curves and specifically examines hardness ra¬ 
tios (Sec. Oi, fractional variability (Sec. 16.2t . and flux-resolved 
spectra (Sec. 1^ . A discussion of the findings (Sec. |7j and final 
concluding remarks (Sec. E complete this work. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 XMM-Newton 

XMM-Newton (Jansen et al. 2001) observations of HE 0436^717 
were serendipitous during the targeted observation of the pulsar 
PSR J0437.4-4711 about 4.2 arcmin away. The first observation 
(hereafter XMMl) was during revolution 0519 on the 9*^ of Oc¬ 
tober 2002 and spanned 70.5 ks. Two MOS detectors (Turner et 
al. 2001) operated in full-frame mode with the thin filter. Due to 
the variable nature of the pulsar target, the pn camera (Striider et 
al. 2001) was set to timing mode also with a thin filter in place. 
The Optical Monitor (OM; Mason et al. 2001) simultaneously col¬ 
lected data from HE 0436-4717 during this observation although 
the object was outside the field of view of the reflection grating 
spectrometers (den Herder et al. 2001; RGSl and RGS2). The OM 
data were limited to observations in the UVWl (Amax = 2675 A) 
and V (Amax = 5235 A) filters. 

A second XMM-Newton observation (hereafter XMM2) was 
during revolution 1835 on the 15*^ of December 2009 and spanned 
130 ks. The EPIC instruments operated in the same manner as 
in XMMl. During XMM2 only the MOS instruments observed 
HE 0436^717 as the source fell outside of the pn field of view 
used for timing mode. 

Data files from both epochs were processed to produce cali¬ 
brated event lists using the XMM- Newton Science Analysis Sys¬ 
tem (SAS) version 12.0.0. The data were examined for background 


flaring and pileup; only minor flaring was seen in the very begin¬ 
ning of XMM2 and those periods were ignored. No pileup was de¬ 
tected as the source was off-axis and the pn was operated in timing 
mode. For the MOS data, source photons were extracted from a 
circular region 35 arcsec in radius and centered on the object. The 
background photons were extracted from an area 50 arcsec in radius 
close to the object and then scaled appropriately. 

Source photons from the pn data were extracted from a rectan¬ 
gular region 2° in width centered on the object and the background 
photons from a larger rectangular region and also scaled. Single¬ 
quadruple events were selected for the MOS data while single and 
double events were selected for the pn; events next to a bad pixel or 
the CCD edge were omitted (i.e. data quality flag set to zero). The 
entire process was then repeated for XMM2. The resulting mean 
count rates were 3.07 count for pn, 1.21 count forMOSl, 
and 1.29 count for MOS2 in the 0.3 - 10.0 keV band. 

Spectra from the MOS data were limited to the 0.5 - 10.0 keV 
range based on cross-calibration uncertainties, while the pn spectra 
were limited to 0.6 - 8.0 keV because of the less accurate calibra¬ 
tion in timing mod^ and high background above 8 keV. 


2.2 Chandra 

The pulsar PSR J0437.4-4711 is a calibration source for the Chan¬ 
dra observatory. Of the eleven observations between 2000 - 2006, 
HE 0436-4717 falls in the field of view during six: ObsIDs 742, 
1850, 6155, 6156, 6157, and 7216. ObsIDs 742 and 1850 were ob¬ 
served with the HRC-I and -S detectors and thus only imaging data 
are available; these were ignored. The ACIS-S data were processed 
using the most recent version of CIAO at the time (v. 4.4.0) and 
preliminary spectra were extracted. After analyzing the spectra it 
was clear that the data were piled up. Several reprocessing meth¬ 
ods were attempted, however the quality of the data were not high 
enough to justify any further analysis. 


^ See XMM-Newton Calibration Technical Note 0018 (XMM-SOC-CAL- 
TN-0018) section 11.2.2 for details 
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2.3 Swift 

HE 0436-4717 has been observed by the Swift X-ray Telescope 
(XRT; Burrows et al. 2005) twice in 2007 and most recently in June 
2013 (hereafter Sw3). The first observation (hereafter Swl) had the 
highest exposure time, followed by the second observation (here¬ 
after Sw2). Processed Burst Alert Telescope (BAT) spectra and the 
diagonal response matrix are taken from the Swift BAT 70-Month 
Hard X-ray Survey (Baumgartner et al. 2014). The BAT data are 
averaged spectra from integration of survey data collected over a 
70-mo period, whereas the XRT observations are individual expo¬ 
sures of a few kiloseconds during pointed observations (see Table 
0 . 

The XRT observations were performed in photon counting 
mode. X-ray data were reduced with the xrtpipeline version 0.12.1 
and spectra were extracted with XSELECT. Optical and UV pho¬ 
tometry was obtained with the UV-Optical Telescope (UVOT; 
Roming et al. 2005). Magnitudes and fluxes were based on the most 
recent UVOT calibration as described in Poole et al. (2008) and 
Breeveld et al. (2010). 

A summary of the data used in this work is shown in Table[T] 


3 PRELIMINARY ANALYSIS OF XMM DATA 

All spectral model fitting was performed using the X-ray spectral 
fitting package XSPEC v. 12.8.1. Model parameters are reported in 
the rest frame of the AON (z = 0.053) and a cosmology of Ho = 
70kms“^ Mpc“^ , qo = 0, and Ao = 0.73 is assumed. All models 
include a Galactic column density of Ah = 1.03x 10^° cm~^ as 
determined from the LAB SurvejQ (Kalberla et al. 2005). Errors 
on model parameters correspond to a 90% confidence level. 

An initial inspection of the XMM-Newton EPIC spectra 
showed consistency between the pn and MOS data for XMMl and 
both MOS spectra for XMM2 within known calibration uncertain¬ 
ties. Since the multi-epoch data — for which only MOS data are 
available — will ultimately be compared, the pn data will not be 
discussed, although it was continually examined for consistency 
throughout the study. The MOS 1 and MOS 2 data at each epoch 
were merged into single spectrum using the FTOOLS program AD- 
DASCASPEC in an effort to increase signal to noise and simplify 
model fitting. 

The merged MOS spectra at each epoch were fit from 2 - 
lOkeV with a power law and the model extrapolated to 0.5 keV. 
The residuals clearly show a soft-excess, common in Seyfert 1 
galaxies, and excess residuals around 6.4 keV (second from top 
panel, Fig.[T](. The addition of a black body at low energies fit the 
soft-excess very well leaving only residuals at higher energies (sec¬ 
ond from bottom. Fig. [TJ. The addition of a narrow (a — 1 eV) 
Gaussian profile at ~ 6.4 keV, common at both epochs, improved 
the fit further but still left significant residuals. Adding a second, 
broader Gaussian profile (E ~ 6 keV and ct « 3 keV) improved 
the quality of fit (Ax^= 14 for 3 additional free parameters; Fig.[T] 
top panel) and residuals (Table[2 Fig.[T] bottom panel). Reversing 
the order of line fitting — adding a broad Gaussian line first — 
does not influence the final fit as the fit statistic remained the same 
(i.e. two Gaussian profiles are still required to reconcile all iron 
residuals). 

The described model provides a good fit to the multi-epoch 
spectra (x^/dof = 1.07/2025; Fig. [T] top and bottom panels). 

^ http://heasarc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 





0.5 1 5 10 


Observed energy (keV) 

Figure 1. Preliminary analysis of the XMM-Newton merged MOS spectra. 
XMM 1 (black) and XMM2 (red) data were fit with a power law plus black 
body and two Gaussians in the 0.5 - 10 keV band (top panel). Extrapolating 
a power law fit to the 2 — 10 keV spectra to lower energies shows clearly 
a soft-excess below about 2 keV and and excess residuals around 6.4 keV 
(second from top). Adding a black body accounts for the soft-excess and 
residuals in the iron band (second panel from bottom). Including a Gaussian 
profile at about 6.4 keV with variable width reduces the residuals further 
(bottom panel). 
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Figure 2. An absorption edge with the same parameter values as those 
found in Wang et al. (1998) analysis of ASCA spectra was added to the 
model and the residuals (data/model) are shown. These data clearly do not 
require the absorption edge and allowing the parameters free to vary re¬ 
turned an absorber with negligible optical depth. XMMl is black, XMM2 
is red. 

Wang et al. (1998) reported an absorption edge around E ~ 
0.83 keV with r ~ 0.26 from their analysis of ASCA spectra. In¬ 
cluding this exact feature in the XMM-Newton spectrum degraded 
the fit (Fig.|2j' Allowing the absorption energy and depth to vary 
did not improve the fit significantly (Ax^= 11 for 2 additional free 
parameters) and the measured edge energy (E = 0.72 ± 0.03) and 
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optical depth (r = 0.06 ± 0.03) were inconsistent with the ASCA 
measurements. Likewise, adding an absorber associated with the 
host galaxy did not enhance the fit and only an upper limit of 
A^h < 10^® cm“^ could be measured. The XMM-Newton spectra 
here do not support significant additional absorption that was re¬ 
ported from the A5C4 analysis (Wang et al. 1998). 

The best-fit phenomological model presented in Table does 
not reveal the physical processes at work in HE 0436-4717. The 
broad Gaussian profile could be associated with blurred reflec¬ 
tion or partial covering absorption. The strength of the soft-excess 
(SE), defined as the flux ratio between the black body and power 
law component between 0.4 — 3keV, could indicate whether the 
low-energy spectrum is sculpted by absorption or reflection (Va¬ 
sudevan et al. 2014), but the ratio measured in HE 0436^717 
(SE « 0.1) is consistent with either process. In the following sec¬ 
tions we examine a physical description for the X-ray emission in 
HE 0436^717. 


4 UV/X-RAY SPECTRAL ENERGY DISTRIBUTION 

X-ray emission provides valuable insight into the accretion process 
in the inner regions closest to the black hole, however UV radiation 
from the accretion disk dominates the bolometric emission. The 
opticalAJV and X-ray bands are investigated together to capture as 
much of the total accretion emission as possible. One objective in 
such broadband analysis is to determine whether or not an AGN is 
potentially bare. 

4.1 The UV-to-X-ray spectral slope 

The slope of a hypothetical power law connecting the X-ray (2 keV) 
and UV (2500 A) spectral bands (ctox, Tananbaum et al. 1979) was 
estimated for HE 0436^717 at both epochs using simultaneous 
XMM-Newton MOS and OM data from the UVWl filter. The UV 
luminosity at 2500 A is ~ 10^® erg . The Qox values were com¬ 
parable at both epochs with Oox.xmmi = -1.27 and aox,XMM 2 = 
-1.26. The higher X-ray flux in XMM2 predicts that the aox slope 
be shallower for that epoch - as is the case. These Qox slopes found 
are reasonable for an AGN of this UV luminosity (Vagnetti et al. 
2013) and all fall within the span of values calculated by Grupe et 
al. 2010 (hereafter GIO) for 92 Type I AGN. The aox values indi¬ 
cate that HE 0436^717 is behaving “normally” (Gallo 2006) and 
extreme absorption or extreme relativistic effects are not likely to 
be present. 

4.2 Detailed SED Modeling Using Swift 

Ideally, one would like to analyze a broadband spectral energy dis¬ 
tribution (SED) that can be used in model fitting directly (e.g. Elvis 
et al. 1994). Swift UVOT data from epoch Swl were used for 
the broadband analysis of HE 0436^717 as it provides data in 
five optical bands (as opposed to XMM-Newton OM, which ob¬ 
served with only two filters). The XRT Swl data provided the 2 
- lOkeVspectrum and the 70-mo BAT data were included for the 
hard X-ray band coverage. 

Following the procedures outlined in Vasudevan & Fabian 
(2009) and then expanded upon by Vasudevan et al. (2009) (here¬ 
after V09), a basic model composed of a broken power law and a 
black body was applied to the SED. The break radius of the bro¬ 
ken power law, Rbr, is defined as the radius at which the photon 
index transitions from its inner value, Ti, to outer, r 2 , in units of 


Rg. Parameters defined in the black body accretion disk model in¬ 
cluded the disk inner radius (Rin), the disk normalization (K), and 
the maximum temperature of the accretion disk (Tmax). Tmax was 
allowed to vary and the inner radius was set to 6 Rg, the innermost 
stable circular orbit (ISCO) of a Schwarzschild black hole. The nor¬ 
malization is defined as X = ' where i is the inclination of 

the disk, Dl the luminosity distance to the source in kpc, and /3 the 
spectral hardening ratio. The black hole mass for HE 0436-4717 
was adopted from GIO as Mbh = 5.9x10^ Mq and a luminosity 
distance of 233.1 Mpc was used0 V09 demonstrated that cos i and 
P could be set to 1 for type-I AGN, and thus the value for K was 
calculated and fixed in the initial fits. 

The assumption regarding is reasonable due to the nature 
of the analysis and the DISKPN model: /3 is the T^oior / T^ffecUve 
ratio that is designed to correct flux in the far-UV band of a spec¬ 
trum in order to preserve the total integrated bolometric luminosity 
(Shimura & Takahara 1995). In the model DISKPN, a decrease in P 
will boost flux across all wavebands, not just in the far-UV. In ad¬ 
dition, as only near-UV data were used in the SED fitting, no color 
correction is ultimately required. As for inclination or cos i, its in¬ 
fluence on the overall normalization is minimal compared to both 
Mbh and Dl, which are both raised to the second power. In gen¬ 
eral, a change in the normalization parameter, K, overall does not 
alter the spectral shape of the component significantly, but rather 
adjusts the peak temperature of the (estimated) black body in order 
to better fit the data. 

Assuming a Comptonization origin for the primary X-rays, 
the low-energy break in the power law corresponded to Tmax 
of the disk component. Galactic absorption was set to A^h = 
1.03x10^° cm“^ and intrinsic extinction was considered in vari¬ 
ous ways, as will be described below. Optical filters (V, B, and U) 
were not used for fitting as host galaxy contamination was not ac¬ 
counted for; the optical data were plotted simply for comparison. 
The XRT data were also fit over a reduced range of 2 - 10 keV. This 
limit was placed to ensure the primary X-ray emission was fit rather 
than the soft excess which may originate from a different compo¬ 
nent. The BAT data were not fit at all and were simply overplotted 
for visual comparison. 

The first model tested (Model A) described the simplest case 
where only Galactic extinction modified the intrinsic broadband 
spectrum. This scenario did not describe the data well (x^ = 102 / 
30 d.o.f), overestimating the UV spectral slope significantly. Galac¬ 
tic extinction alone can not explain these data. Host galaxy extinc¬ 
tion was then added to the model and allowed to vary (Model B). 
The fit improved by Ax^= 71 for an additional free parameter. The 
extinction that described the SED best was Eb-v = 0.095±0.01L 
This scenario describes the broadband spectra very well, however it 
seems contradictory with the X-ray analysis performed previously. 
Assuming the estimated extinction, the corresponding column den¬ 
sity can be calculated using the relationship outlined by Predehl & 
Schmidt (1995): Vh = 5.3 x 10^^ cm“^£B-v. For £b-v rang¬ 
ing from 0.084 - 0.105, A^h = (5.04±0.58)x 10^° cm-^ which is 
at least 50 x the upper limit measured from the X-ray spectra (see 
Section[3. Moreover, GIO measures extinction from the Baimer 
decrement and finds Eb-v = 0.04, which is much lower than the 
intrinsic extinction estimated in Model B. 

Next, the intrinsic extinction measured from the Balmer decre¬ 
ment was adopted and the K was allowed to vary (Model C). 
The best value was K = (52.16±14.30) xl0“"' M 0 ^kpc“^ and 

® http://www.astro.ucla.edu/ wright/CosmoCalc.html 
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Table 2. The best phenomenological model fit to the merged XMM-Newton MOS spectra. The model components and model parameters are listed in Columns 
1, and 2, respectively. Columns 3 and 4 list the parameter values during the XMM1 epoch and XMM2 epoch, respectively. Fluxes are corrected for Galactic 
absorption and are reported in units of ergcm“^ . The soft-excess flux is measured between 0.4 — 3.0keV. Parameters linked across the epochs are 
denoted by dots and those without en'ors are frozen at the given values. 


(1) 

Model Component 

(2) 

Model Parameter 

(3) 

XMMl 

(4) 

XMM2 

Black body 

fcr(eV) 

F'oA-akeV 

77 ±6 

0.68±0.05 X 10-*2 

79 ±5 

0.72±0.05 X 10-12 

Power law 

r 

P'oA-akeV 

2.14 ±0.04 
5.36±0.05 X 10-12 

2.23 ±0.03 
6.71±0.05 X 10-12 

Narrow 

Gaussian 

E{keV) 
a (keV) 

EW (eV) 

6.42 ± 0.06 

1 

46 

44 

Broad 

Gaussian 

T;(keV) 
a (keV) 

EW (keV) 

^•^^-3.36 

2.5 

2.3 

Observed 

Flux 

F'2-iokeV 

P'o.5-10keV 

4.30 X 10-12 

8.22 X 10-12 

4.65 X 10-12 

9.52 X 10-12 

Soft- 

Excess 

Tbb ! T’pl 

0.13 ±0.01 

0.11 ±0.01 

Fit Quality 

xS/dof 

1.07/2025 



Table 3. The four models used to fit the UV-to-X-ray SED. Column 1 lists the model parameters and goodness-of-fit. Column 2 shows the parameter values 
for each model. Parameters that are fixed are listed without errors. 


(1) 

Parameters 

A 

(2) 

SED Models 

B C 

D 

001-100 keV 

6.39 

15.17 

6.47 

4.99 

(Xl0‘*‘*ergs 1) 


0 

0.095±0.011 

0.04 

0 

EB-V,Gal 

0.009 

0.009 

0.009 

0.009 

Tmax (eV) 

5.26±0.01 

6.85±0.21 

3.11±0.25 

2.35±0.14 

X(xlO'*M02kpc-2) 

6.41 

6.41 

52.87±14.30 

107.54±25.77 

power law norm. (xl0~^) 

1.36±0.06 

0.63±0.06 

6.22±1.42 

14.16±2.59 

(ph ke'V 1 cm 2 s 1 ) 

TboI ! TBdd 

0.08 

0.20 

0.09 

0.07 

xV d.o.f. 

3.50 / 30 

1.06/29 

1.03/29 

1.02/29 


^ 0 . 001-100 kev = 6.47x10'*'* ergs“*, giving an Eddington lumi¬ 
nosity ratio of LboI / ^Bdd= 0.09. Model C provides a better fit 
than Models A and B, and the treatment of the extinction is physi¬ 
cally motivated by the measured Balmer decrement. 

The final possibility tested is that the SED can be entirely de¬ 
scribed by allowing only the disk normalization to vary with only 
Galactic extinction (Model D). The fit statistic was the compara¬ 
ble to Model C and a: = (107.54±25.77)x 10* Mo^kpc”^ In 
this case, Lq.ooi-ioo kev = 4.99x 10*“* erg s“* and so LboI / iEdd= 
0.07. 

A summary of the model variations with their respective pa¬ 
rameters and Eddington ratios are listed in Tabled The SED mod¬ 
els with high extinction values seem at odds with the initial X-ray 
analysis, extinction measured from the Balmer decrement (GIO), 
and the UV spectral analysis of Halpern & Marshall (1996). Alter¬ 


natively, the SED can also be explained if we assume the source 
is not substantially reddened, but that the ratio (Mbh / is 

uncertain hy a factor of 6 - 10. This is perhaps not unreasonable. 
The black hole mass for HE 0436^717 was not from reverberation 
mapping, but rather based on a single-epoch mass estimate which 
could be uncertain by a factor of a few (Vestergaard & Peterson 
2006). 

The Swift BAT data were included for visual inspection of the 
high energy trend, however the BAT data were not fitted. Although 
the uncertainty in the BAT data allow for a wide spread, the ex¬ 
trapolated models were all consistent with these data. Model C (i.e. 
Balmer decrement based intrinsic extinction and K allowed to vary) 
is shown for both the observer and source frame (see Fig.^. 

The Eddington ratios determined from the Swift SED analysis 
agree with the value of 0.08 found by GIO for this object. An Ed- 
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Figure 3. An SED comparison of Model C with intrinsic extinction de¬ 
termined by the B aimer decrement and the disk parameter K allowed to 
vary. The solid lines represent the total model in the source (red) and ob¬ 
served (blue) reference frames. The dashed line represents the disk com¬ 
ponent while dotted shows the primary X-ray component. Only the XRT 
data (black squares) from 2.0 - lO.OkeV and the UV data of the UVOT 
(tiny maroon circles) were fit, the 70-month Swift BAT data (green trian¬ 
gles) were overplotted for visual inspection. Only the XRT data that were 
fit (i.e. between 2-10 keV) were included for clarity. 

dington ratio around 0.08 is not unusual, although it does fall on 
the lower end of the probability distribution described by V09 in 
their analysis of 26 low-z AGN (see V09 Fig. 15, lefthand panel). 
A lower Eddington ratio implies HE 0436-4717 is not a rapid ac- 
creter. 


5 DETAILED X-RAY SPECTRAL ANALYSIS 

The first glance at the broadband X-ray spectra of HE 0436-4717 
showed the presence of multiple components and the nature of the 
excess residuals above ~ 6 keV remains unclear. Three scenarios 
are considered to describe the X-ray spectra of HE 0436^717 in a 
self-consistent manner: partial covering absorption, blurred reflec¬ 
tion, and soft Comptonization. 

5.1 Partial Covering Absorption 

Turner & Miller (2009) explain in detail how gas with a range of 
column densities can influence the primary component of an AGN. 
As absorbers move in the field of view, some of the intrinsic radi¬ 
ation is both scattered and absorbed. Both changes in the covering 
fraction and column density can create spectral variability over a 
range of timescales. A strong soft excess can appear at lower en¬ 
ergies and even the Fe Ka line can be broadened by absorption 
with the red shoulder being enhanced by attenuation rather than by 
relativistic effects. A thorough attempt was made to fit the X-ray 
spectrum of HE 0436-4717 with absorption in an effort to deter¬ 
mine the origin of the X-ray emission in this object. 

The model ZPCFABS was used to begin testing for possible 
neutral absorption signatures. A single, neutral absorber was ap¬ 
plied to a simple power law and fit to both epochs of the merged 
MOS data. Scenarios were tested in which only the covering frac¬ 
tion remained free to vary between the epochs (i.e. variability at¬ 
tributed to absorbers) and also having only the primary compo¬ 


nent vary instead (i.e. variability attributed to primary emitter). The 
overall goodness-of-fit remained poor in both cases. Double neu¬ 
tral absorption was an improvement, though the model continued 
to underestimate the data at higher energies. This scenario called 
for the two absorbers to have very different densities (see Table 
0 and could potentially describe a single absorber with a density 
gradient along the line of sight, rather than two completely inde¬ 
pendent bodies. The column densities remain consistent in the two 
epochs and so would not alone explain any spectral variability seen 
over longer timescales. Allowing only the primary emitter to vary 
instead did not improve the fit quality significantly. Apparently, if 
there is any variability in the primary emitter on long timescales, it 
is subtle. 

5.2 Ionized Partial Covering 

Next, a fit with a single ionized absorber was attempted. The ion¬ 
ized absorber was modeled by ZXIPCF. Once again, all parameters 
except for covering fraction were linked at the epochs in an effort 
to describe the long-term variability with changes in an absorber. 
There were residuals to the model and the fit quality was poor (x^ 
/ d.o.f. = 2415 / 2032). The covering fraction was not significantly 
different from one epoch to the next. 

If the photon index of the primary component is allowed to 
vary while the absorber is frozen and linked across the epochs, the 
fit improves by Ax^ = 213 for two additional free parameters, how¬ 
ever some positive residual features remain. The photon index in 
XMM2 increases, going from F = 2.16±0.01 to F = 2.24±0.01, 
consistent with the trend found in the toy model fits. It appears that 
the model with a single ionized absorber detects a slight, long-term 
variability from the primary emitter while the model with a single 
neutral absorber does not. 

A second ionized absorber was added to the model and the 
parameters reset so that only the covering fractions were unlinked 
across the epochs. The fit was significantly better than that of the 
single absorber (x^/ d.o.f. = 2150 / 2028) and no residual features 
remained (see Fig. 0, although the covering fractions for both ab¬ 
sorbing bodies remained similar at both epochs (see Table0. Thus, 
as a final exercise, one absorber was kept constant throughout the 
epochs while the other was allowed to vary in Ah and with F 
free to vary as well. The goodness-of-fit remained unchanged and 
the parameters that were allowed to be independent in XMM2 did 
not significantly change from their XMMl values. 

None of the ionized absorption model variations defined the 
source of long-term variability uniquely, however, the scenario in 
which only the covering fraction varied between epochs did result 
in the best goodness-of-fit. The most satisfying absorption model 
overall describes two ionized absorbers, the first more ionized and 
with a higher column density than the second, with both having 
comparable covering fractions. Both absorbers are found to have 
redshifts comparable to the source that remained steady across the 
epochs. 

The ionized partial covering model was extrapolated to higher 
energies and visually compared to the Swift BAT data. Model resid¬ 
uals show that, despite the partial covering scenario’s ability to de¬ 
scribe the 0.5 - lOkeVspectra well, it under-predicts the X-ray flux 
above lOkeV (Fig.0 residuals). However, photoelectric absorption 
will also produce its own reprocessed emission in addition to the 
intrinsic radiation it attenuates (e.g. Turner & Miller 2009). Repro¬ 
cessed emission from absorption would affect both the hard X-ray 
flux and flux of the Fe Ka line. 

Following Reynolds et al. (2004) we estimate the strength of 
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Table 4. The best spectral absorption models in the 0.5 - 10.0 keV band. Models are listed in Column 1 and their components in Column 2. Model parameters 
are in Column 3 and their values in Columns 4 and 5 for XMMl and XMM2, respectively. Parameters linked across the epochs are denoted by dots. 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

Model 

Model Component 

Model Parameter 

XMMl 

XMM2 

double neutral 

absorber 1 

Nh (cm- 2 ) 

26.30+? ^xl022 


absolution 


CF 

0.61 ± 0.02 

0.56±0.02 


absorber 2 

Nh (cm -2 ) 

1.33+°i°x1022 




CF 

0.48±0.02 

0.42±0.02 


primary flux 

r 

2.66±0.03 



fit quality 

xj/d.o.f. 

1.11/2030 


double ionized 

absorber 1 

Nh (cm- 2 ) 



absorption 


log( 0 [erg cms“l ] 

1 

±.oy_Q 3Q 




CF 

0 57+0 02 
-0.03 

0.63±0.03 


absorber 2 

Nh (cm- 2 ) 

1.75t° ^ 3 x 1022 




log (0 [ srg cm s“ ^ ] 

0.48t°;« 




CF 

0.52±0.03 

0.44±0.04 


primary flux 

r 

2.63±0.05 



fit quality 

Xg/d.o.f. 

1.06/2028 




Observed energy (keV) 

Figure 4. The ionized partial covering model (Table |4j and its relation to 
the MOS spectra. A pair of ionized absorbers were fit from 0.5 - 10.0 keV 
with the covering fraction as the only parameter varying at both epochs. The 
XMMl (black) and XMM2 (red) data are shown with the constant intrinsic 
power law determined by the model (dotted blue) in the top panel. The 
model was extrapolated above 10 keV for visual comparison with 70-month 
Swift BAT data (yellow) and the corresponding residuals are shown (bottom 
panel). 


an Fe Ka emission line that would be produced from the denser 
of the two absorbers in the partial covering model. Assuming the 
number of ionising photons between 7-20 keV that are reemitted 
in the neutral iron line at 6.4 keV is governed by the fluorescent 
yield of neutral iron (0.347) we estimate the flux in the iron line 
should be approximately 2.70x10“'^® ergcm“^ for XMMl, 
2.55 X10“^® erg cm“^ for XMM2. Such a narrow line would 
have an equivalent width of EW « 210 eV and 200 eV, in the 
spectrum of XMMl and XMM2, respectively. The strength of the 
narrow Fe Ka emission line in the spectra of FIE 0436^717 is only 
about 40 eV, much weaker than expected if arising from reprocess¬ 


ing, but we note that various effect such as Compton-thick lines-of- 
sight could generated weaker lines than prediced (e.g. Yaqoob et al. 
2010 ). 

Finally, an unblurred reflection component REFLIONX (Ross 
& Fabian 2005) was added to the absorption model to model for 
any narrow Fe Ka emission, but also to adjust the residuals above 
10 keV seen in Fig.|4] The modification above 10 keV was negligi¬ 
ble, which is expected given the steep photon index required to fit 
the spectra below 10 keV. 

5.3 Blurred Reflection 

The multi-epoch observations of HE 0436^717 were next fit with 
a blurred reflection model. The primary X-ray source was modelled 
as a cutoff power law with the high-energy cutoff fixed at 300 keV. 
The backscattered emission from the accretion disk, modelled with 
REFLIONX, was modified by relativistic effects close to the black 
hole (Fabian et al. 1989). We adopted KDBLUR2 for the blurring 
kernel, which incorporates a twice broken power law for the emis- 
sivity profile. 

The disk inclination (i), iron abundance (Afb), and inner disk 
radius (Rin) are not expected to change over the observable time 
scales and were linked between epochs. The iron abundance was 
initially fixed to solar, the outer disk emissivity profile (qout) was 
fixed to 3, and the break radius (Rb) where the emissivity index 
changes from qi„ to qout was fixed to Rb =6 Rg = 6 ^4^. At each 
epoch, the ionization parameter = ‘iitFjn where n is the hydro¬ 
gen number density and F is the incident flux), normalization of the 
reflector, the incident power law photon index, and normalization 
were all free parameters. 

The measurement of iron abundance is rather sensitive and de¬ 
generacies have been found with other fit parameters like emissiv¬ 
ity, inclination, and black hole spin (or comparably, the inner disk 
radius). Contour plots were generated between various fit parame¬ 
ters to determine if the low iron abundance could be an artifact of 
the modelling, but in all cases the iron abundance remained consis¬ 
tently low (Fig. Is}. 

A second REELIONX component that was not blurred and 
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Table 5. The best-fit blurred reflection model fit simultaneously to the XMM-Newton merged MOS spectra. The model components and model parameters are 
listed in Columns 1, and 2, respectively. Columns 3 and 4 list the parameter values during the XMMl epoch and XMM2 epoch, respectively. Parameters that 
are fixed ai‘e listed without en'ors and the dots denote those that are linked across the epochs. The reflection fraction (7^) is approximated as the ratio of the 
reflected flux over the power law flux in the 0.1-100keV band. Fluxes are corrected for Galactic absorption and are reported in units of ergcm”^ . 


( 1 ) 

(2) 

(3) 

(4) 

Model Component 

Model Parameter 

XMMl 

XMM2 

Incident 

r 

2.03±0.02 


Continuum 

Bout (keV) 

300 



Gj. 1-100 keV 

1.35±0.06 X 10-11 

1.66±0.08 X 10-11 

Blurring 

9in 

« q+2.0 

0-0_Q 7 

5.6 ±0.7 


^in ^g) 

< 1.80 



-Rout ( -Rg ) 

400 



Rb ( Rg ) 

6 



Qout 

3 



in 

45l« 


Reflection 

^ (ergcms”^ ) 

150t“0 

56l« 


Ape (solar) 

0.36t°°^ 



Rq. 1-100 keV 

4.90±0.63 X 10-12 

4.79±0.52 X 10-12 


n 

0.36±0.05 

0.29 ±0.03 

Distant Reflector 

^ (ergcms”^ ) 

1.0 



Ape (solar) 

1.0 



r 

1.9 



Gl.1-100 keV 

1.621“;“ xlO -12 


Fit Quality 

x 2 /dof 

1.06/2024 






Inner emlssivity index 


Figure 5. Contours are calculated to examine the dependency of the measured iron abundance with various other parameters. The iron abundance is compared 
to the inner disk radius (left), disk inclination (center), and inner emissivity profile (right). In the right panel the solid and dashed contours are for data set 1 
and 2, respectively. The contours correspond to a 2.3, 4.61 and 9.21 delta fit statistic. 


whose flux remains linked between epochs was also included 
to mimic emission from a distant reflector like the torus. The 
ionization parameter of this distant reflector was fixed at ^ = 
lergcms”^ and the iron abundance set to the solar value. The 
photon index of the power law source illuminating the torus is fixed 
at r = 1.9, the canonical value for AGN (e.g. Nandra & Pounds 
1994). There is no reason to assume that the radiation incident on 
the distant reflector will be the same as that from the primary com¬ 
ponent and if a value of T = 2.0 is used instead it does not signifi¬ 
cantly alter the fit statistic. 

The described model produced a reasonable fit to both spectra 
(X?/dof =1.08/2025), but left excess residuals above 5keV. Per¬ 
mitting Qout and Rh to vary at each epoch did not substantially 
improve the fits, but a significant improvement (Ax^=33 for 1 ad¬ 
ditional fit parameter) was found if the iron abundance was allowed 
to vary. The resulting fit was good (Xy/dof=l.06/2024; Table O, 
but noteworthy was the low, sub-solar iron abundance measurement 
0.361°°^). 

Further examination to determine if the low iron abundance 


was a result of the fit procedure were carried out. The iron abun¬ 
dance of the distant reflector was linked with the inner, blurred re¬ 
flector to determine if enhancing the iron “signal” could result in a 
more typical measurement. This fit was an improvement over the 
initial fit, with the abundance fixed at solar (Ax^=25 for 1 addi¬ 
tional fit parameter), but the resulting measurement was compara¬ 
ble to the measurement from the blurred component alone {Ape = 

0.33t°;i^). 

Strong degeneracy has been noted between the black hole spin 
parameter (a) and the iron abundance (e.g. Reynolds 2013). The 
weak reflection features and moderate signal-to-noise of the data do 
not justify a thorough examination into black hole spin. However, 
the possibility that the low iron abundance may be rectified with 
treatment of the black hole spin should be considered. With this in 
mind, the blurring kernel KDBLUR2 was replaced with RELCONV 
(Dauser et al. 2010) which includes (retrograde and prograde) black 
hole spin as a measured parameter. The fit quality was comparable 
to that achieved with KDBLUR2 (x^/dof=l. 06/2024). The spin pa¬ 
rameter was high (a « 0.96) and was consistent with the value 
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Observed energy (keV) 

Figure 6. The bluiTed reflection model for the two epochs (top panel). The 
primary X-ray component is dotted while the bluiTed reflector is dashed. 
The distant reflector (same in both epochs) is solid blue. The total summed 
model is solid black and red, for XMMl and XMM2 respectively. The 
model was extrapolated above lOkeV for visual inspection with the BAT 
data (yellow) and the residuals (model/data) are shown (bottom panel). 

of the inner disk radius reported in Table |5] The iron abundance 
continued to remain low (Aps ~ 0.28). Forcing the model to ac¬ 
commodate a retrograde spin (—0.998 ^ a ^ 0) only resulted in 
a poorer fit (xJ/dof=l.l 1/2024) and did little to enhance the iron 
abundance (Ape « 0.56). 

Further consideration was to adopt a different reflection code. 
The reflection model REFLIONX was replaced with XILLVER (Gar¬ 
cia et al. 2013). There are several differences between XILLVER and 
REFLIONX including XILLVER’s ability to calculate reflected inten¬ 
sity for each photon energy, position, and viewing angle whereas 
REFLIONX uses an angle-average integration for reflected intensity. 
Primarily, we consider if these differences could be attributed to 
the low iron abundnace in HE 0436-4717. Substituting REFLIONX 
for XILLVER in the original model (Table|5]( resulted in an equally 
good fit (x^/dof=l.07/2024) with comparable parameters mea¬ 
surements. We note the iron abundance boundary value in XIL¬ 
LVER is 0.5, which makes it impossible to explore values as low 
as we found with REFLIONX. Nonetheless, the best-fit iron abun¬ 
dance pegs at the lower limit of the model (4_Fe=0.50(^Q ^®). 

5.4 Soft Comptonization 

Both the absorption and reflection scenarios addressed in the previ¬ 
ous sections fit the soft-excess of HE 0436^717 well and it is sta¬ 
tistically impossible to distinguish between the two scenarios with¬ 
out further analysis. A third possibility is that the soft-excess could 
be emerging from a second Comptonizing region. 

Observations of binary black holes have shown that spectral 
fitting can require two Comptonizing regions (Kubota & Done 
2004 and Done et al. 2012, hereafter D12). While there are sev¬ 
eral possibilities for the source of the second softer Compton up- 
scattered photons, D12 assume such soft Comptonization is origi¬ 
nating from the accretion disk itself. The authors postulate that disk 
structures could deviate from the standard Shakura-Sunyaev model 
at inner radii, as has been proposed for some of the more extreme 
Seyfert galaxies like Narrow-Line Seyfert Is (NLSls). In particular 



0.5 1 5 10 too 

Observed energy (keV) 


Figure 7. The double-Comptonization model for the two epochs (top 
panel). The primary X-ray component is dotted while the second Comp¬ 
ton component is dashed. The distant reflector (same in both epochs) is 
solid blue. The total summed model is solid black and red, for XMMl and 
XMM2 respectively. The model was extrapolated above lOkeV for visual 
inspection with the BAT data (yellow) and the residuals (model/data) are 
shown (bottom panel). 

they consider enhanced dissipation of disk outer layers that would 
create a Compton up-scattered flux in the soft X-ray band. 

As an initial investigation into a double-Comptonization sce¬ 
nario, the merged MOS spectra were fit from 0.5 - 10.0 keV with a 
basic power law -l- Comptonization model with the second Comp¬ 
tonization component provided by COMPTT (Titarchuk 1994). The 
seed photons were assumed to be supplied by the intrinsic UV disk 
radiation below the observed bandpass and both the seed photon 
temperature and the temperature of the second Comptonizing re¬ 
gion were allowed to vary. This lone double-Comptonization model 
did not fit well (x^ / d.o.f. = 2319 / 2030) and there were clear 
residuals around 6.4 keV. When a distant reflector was added, the 
fit improved by Ax^ = 99 for two additional free parameters, but 
the new model did not account for all of the iron line flux. If a broad 
line was added in the form of blurred reflection, the spectra were 
well-modeled in the 5 - 8keV range, however the normalization 
of the second Compton component dropped by an order of magni¬ 
tude and the optical depth reduced from around 2 keV to r ~ 0.09 
as blurred reflection now accounted for the soft band. A simple 
double-Comptonization scenario appears to model the soft-excess 
appropriately, but does not fit some of the curvature between 0.5 - 
10.0 keV. 

In an attempt to model the full 0.5 - 10.0 keV spectra more 
physically, the data were fit with OPTXAGNF (D12). Black hole 
mass and luminosity distance were frozen at their known values 
of 5.9x10^ Mq and 233.1 Mpc, respectively. The Eddington frac¬ 
tion was set at log(LILEdd) = -1.05 as determined by the SED fit¬ 
ting, a value that — as mentioned previously — agrees with the 
fraction found by GIO. Black hole spin was set to a = 0.5 to start 
and the normalization was fixed at unity as per model instruction^- 
For the initial fitting, only the power law parameters were allowed 
to vary independently across the epochs and all other parameters 
were frozen. 

https://heasarc.gsfc.nasa.gov/xanadu/xspec/models/optxagn.html 
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Table 6. Summary for the best soft Comptonization model. All parameters 
are frozen with the exception of: Rear, kTe, t, T,fpi and all parameters are 
linked across the epochs with the exception of F. 


(1) 

(2) 

(3) 

(4) 

Component 

Parameters 

XMMl 

XMM2 

Distant Reflection 

Ape (solar) 

1.0 



Fdr 

1.0 



^(ergcms“^ ) 

1.0 



z 

0.053 



norrYif^^ 

(1.56 ±0.35) X 10-^ 


Absorption Edge 

Badge (keV) 

ft 

‘ ' -0.02 



’^max,edge 

(9.13 ±0.02) X 10-^ 


Soft Comptonization 

Mbh (M0) 

5.9x10’’ 



Dl (Mpe) 

233.1 



log(L/ Lpdd) 

-1.05 



logi?out 

3 



a 

0.8 



^cor ( H-g) 

-|-| t:rc + 5.41 



kTe (keV) 

0 g1+0-20 

0.11 



T 

7 70 + 1.19 
' •''^- 1.27 



r 

-l--oU_o_o2 

-•-•yo_o.o5 


/pi 

0.44±0.01 


xV d.o.f. 

1.09/2029 




Parameters were systematically thawed and unlinked to grad¬ 
ually increase complexity as needed to improve the fit. The param¬ 
eters /pi, 7?cor, kTe, and remained consistent between the epochs 
when allowed to vary and so were kept linked. Although it is rea¬ 
sonable to allow these components to vary between epochs, they 
did not improve the quality of the fit if unlinked and so, after an 
initial test of allowing them free to vary at both epochs, they were 
re-linked. 

The initial fit was very poor (x^ / d.o.f. = 2475 / 2032) with 
a clear iron line residual and other residual features at about E ^ 
1 keV. The fit proved difficult to constrain when the spin parame¬ 
ter was allowed to vary freely, so the fit quality was instead tested 
at different fixed values of a. A zero spin only worsened the fit 
(Ax^ = 148), while testing a maximum spin value of a = 0.998 
improved the fit by Ax^ = 303 for no additional free parameters. 
Fixing the spin at a moderately high value of a = 0.8 improved the 
goodness-of-fit even further, with Ax^ = 12 for the same number of 
free parameters, and a was kept frozen at this value. If the Edding¬ 
ton ratio is allowed to vary with spin, neither parameter changes 
significantly from their best-fit values and the fit statistic remains 
comparable. 

To improve the fit, a distant reflector was added in the same 
manner as Section [53] The distant reflector improved the fit signif¬ 
icantly, Ax^ = 101 for one additional free parameter, and the iron 
line residuals were resolved. Curiously, the data remained overes¬ 
timated around 0.8 keV until an absorption edge was included at 
the redshift of the source, with both its energy and optical depth 
free to vary. The absorption feature was found to have an energy 
of £ ~ 0.77 keV and r 0.09 - consistent with O VII, though the 
accompanying O VIII absorption edge was not detected. We noted 
this is the only model that shows marked improvement with the ad¬ 
dition of an absorption edge. The addition of both a distant reflector 
and an absorption edge provided the best description of the data in 
the soft Comptonization scenario (x^ / d.o.f. = 2211 / 2029) and 
the model that most satisfies the data is summarized in Table|^ 

The model was extrapolated into the 10 - 100 keV band and 
the BAT data was over plotted for comparison. The model compo¬ 
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Figure 8. The soft Comptonization model (Table is shown in the top 
panel. XMMl is black and XMM2 is red. The hard power law component 
is shown as the dashed line. The soft Comptonization component that is 
common at both epochs is shown as a dotted blue curve. The distant re¬ 
flect, also common at both epochs, is shown as the fine, solid, blue curve. 
The model was extrapolated >10keV for comparison with the BAT data 
(yellow) and the con'esponding residuals (data/model) are shown (bottom 
panel). 


nents (top panel) and residuals that include the BAT data (bottom 
panel) are shown in Fig. [S] The BAT data were only over plotted 
for comparison and were not included in the fit. The model slightly 
underestimates the fluxes seen above lOkeV. 

Unlike the partial covering or blurred reflection models, the 
OPTXAGNF model provides some information regarding the UV, 
disk dominated emission. Thus, it is advantageous to explore if the 
best soft Comptonization model is consistent with the Swift UVOT 
data. Simply extrapolating the model to the UV band shows a poor 
fit using the assumed black hole mass and luminosity distance, as 
did our initial SED fits (Section [4.2t . Allow the black hole mass 
and distance to vary does not significantly improve the fit since 
the statistics are driven by the X-ray data, but it does suggest there 
could be a discrepancy of at least 2 in the Mbh/77l ratio compara¬ 
ble with our finding with the SED modelling. 

In an further attempt to fit the UV data, the parameters defin¬ 
ing the disk component — black hole mass, Eddington ratio, and 
spin — were all allowed to vary. The fit is comparable to that of the 
X-ray only data (x^ / d.o.f. = 2246 / 2029), however the model un¬ 
derestimated the UV flux. The black hole mass increased by ~ 27% 
while spin decreased by ~ 20%; the Eddington ration remained un¬ 
changed. If the remaining parameters - those that describe the X- 
ray components and the absorption edge - are free, the fit statistic 
does not change significantly and neither do the parameter values. 
Lastly, an absorber was added to the model to mimic intrinsic ex¬ 
tinction in the UV. This broadband soft Comptonization model with 
absorption fit well (x^ / d.o.f. = 2213 / 2028), however it predicted 
an intrinsic extinction of EB-v,int = 0.055 which requires a col¬ 
umn density of Ah = 2.99 x 10^° cm“^ : 29x higher than that pre¬ 
dicted by the X-ray data (Ah< 10^® cm“^ ). Once again it appears 
that an absorption scenario is statistically acceptable for HE 0436- 
4717, albeit inconsistent with further scrutiny. 
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Bin time: 500.0 s 



Bin time: 500.0 s 



Figure 9. MOS light curves for both XMMl (left) and XMM2 (right) binned 
displays more variability. Zero marks the start of each exposure. 

6 TIMING ANALYSIS 

Source photons were extracted to create light curves in a variety of 
energy hands between 0.2 - 12.0 keV, including a broadband 0.2 - 
12.0 keV bin. The light curves from all EPIC instruments that col¬ 
lected data proved consistent in the XMMl and XMM2. All subse¬ 
quent timing analyses utilized MOS data only. MOS 1 and 2 light 
curves were combined at each epoch, taking care to match the ob¬ 
servation start and stop times. The 0.2 - 12.0 keV merged MOS 
light curve for each epoch is shown in Fig.|3 


6.1 Hardness Ratios 

Ratios of softer energy band count rates over those of a harder en¬ 
ergy band were calculated for four energy band pairs between 0.2 
- 12.0keV range. Bins were selected such that there were a com¬ 
parable number of counts in each band. These hardness ratios (HR) 
were analyzed two ways at each epoch: flux-dependent HR; and as 
time-dependent, normalized HR. 

The hardness as a function of count rate did not reveal any sig¬ 
nificant correlation at either epoch. This does not mean the spectra 
are not variable, but only that any variability is not significantly de¬ 
pendent on the specific flux state. The hardness as a function of time 
was similar: within the calculated uncertainties it appears there is 
little time-dependent variation in the soft bands over time. All HR 
trends for the X-ray data were nearly linear within uncertainties and 
so indicate little to no change in hardness ratio with respect to ei¬ 
ther flux or time. This does not mean that the spectral components 
are not variable, but rather that any variability that may be present 
occurs in such a way as to maintain a constant flux ratio between 
energy bins. A more thorough analysis is needed to further probe 
the temporal behavior in HE 0436^717 and this is explained in the 
following sections. 


6.2 Fractional Variability 

The fractional variability amplitude (Edelson et al. 2002), TVar, is 
a useful measurement used to quantify the intrinsic variability of 


at 500 s each. XMM2 had nearly double the duration (note x-axes) and also 


a light curve while also accounting for uncertainty. This parame¬ 
ter compares the amplitude of any variations present in each en¬ 
ergy band to determine if the amplitudes in one energy band are 
larger than those in another. The F^ai thus compares the relative 
“strengths” of observed variability across two energy bands, inde¬ 
pendent of time or flux. The i^var is defined as, 

i^var = (alO (1) 

where (X) is the mean count rate, the total variance of the light 
curve, and (aerr) the mean error squared. Uncertainties were deter¬ 
mined following the procedure of Ponti et al. (2004). 

The TVar was calculated for both epochs using merged MOS 
light curves. The light curves were binned at 1000 s and a total of 
nine energy bins were used spanning 0.2 - 10.0 keV; data above 
lOkeV were ignored due to the low signal-to-noise. XMM2 dis¬ 
plays a higher TVar than XMMl overall, which is apparent in the 
light curves (e.g. Eig.|^. 

In terms of the blurred reflection model, variations in the pri¬ 
mary component (i.e. power law) may solely account for the shape 
of the Fva,T spectrum. Erom the blurred reflection model found pre¬ 
viously, the fraction that the power law contributes to each energy 
band was determined and overplotted on the TVar spectrum for 
both epochs (Eig. [TO] dashed blue line). In the blurred reflection 
model it seems reasonable to attribute most of the rapid variabil¬ 
ity to changes in the power law normalization, specifically. In this 
very simplified scenario, the reflection model prediction disagrees 
with the data in a few bins (around 0.4keVand between 1 - 2keV, 
for example). The best that can be said is that the overall model 
trend does agree with the linear shape of the fractional variability 
spectrum. 

When the process is repeated for the partial covering absorp¬ 
tion model, the predicted FVar profile is less constrained (Eig. [TO] 
dotted green line). Testing the partial covering scenario is more 
challenging as the source of variability could be the absorber as 
well as the primary component. The same difficulty lies in mod¬ 
eling soft Comptonization: there could be many sources of vari¬ 
ability. Eor completeness, a profile was created using the double- 
Comptonization model of Section lS^ with the hard Compton com- 
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Figure 10. Fractional variability analysis with spectral model predictions. 
The fraction that the power law contributes to each energy band is overplot¬ 
ted on the FVar spectrum for both XMMl (top panel) and XMM2 (bottom 
panel). The reflection (blue dashed line), absorption (green dotted line), and 
double-Comptonization (magenta dot-dashed line) model predictions are 
compared for each epoch. 

ponent assumed to be the sole source of variability (Fig. [TOl dot- 
dashed magenta line). 

Looking at the three model predictions for both epochs, it is 
clear that there must be more to the story. Perhaps only a fraction of 
the total power law flux is varying or perhaps an additional compo¬ 
nent contributes to the total variability. The TVar analysis indicates 
that multiple parameters need to vary if any model is to describe 
the spectral data. 

6.3 Flux-Resolved Spectra 

Flux-resolved spectra were made in an effort to resolve the mecha- 
nism(s) behind short-term variability in HE 0436^717. As XMM2 
was considerably longer and shows larger amplitude variations, 
data from the second epoch alone were used to create X-ray spectra 
with < 3 count and > 3 count . 

A difference spectrum (high flux - low flux) was created to 
identify components that vary between the flux states. If changes 
between the flux levels were only due to changes in the normal¬ 
ization of the spectral components then a single power law should 
describe the difference spectrum well. Fitting the 3 - lOkeV band 
with a power law (F = 2.00 ± 0.16) resulted in a good fit. Ex¬ 
trapolating the power law to 0.5 keV reveal a strong excess toward 
lower energies. Normalization changes may be responsible for the 


majority of variability in the broadband spectra, although changes 
in the shape of the spectral components at lower energies may be 
important. 

The two flux-resolved spectra were then fit with the models 
for partial covering absorption, blurred reflection, and Comptoniza- 
tion to examine the origin of the rapid variability. In the absorption 
model, changes in X-ray flux could be caused by either the pri¬ 
mary emitter or by the absorbing material. Two sources of ionized 
absorption explain the X-ray spectra best (Section ED. one with 
higher column density than the other. Geometrically, these could 
represent two distinct absorbing materials or a single absorber with 
a density gradient. Ionization could also vary across a single larger 
body as the material closer to incident radiation would experience 
more illumination than that on the far side. As absorbing material 
moves, overlaps, or condenses/disperses a spectrum can change. 
The simplest scenario assumes variability caused only by a change 
in covering fraction, with all other properties of the absorbers re¬ 
maining constant. 

Allowing only the covering fractions to vary between flux 
states in the ionized absorption model from Section IF!2l returned 
d.o.f. = 909 / 768 and moderate residual features. There was 
a small discrepancy in the spectra around 5 keV where the model 
over-predicted the low flux state and under-predicted the high flux 
state. In addition, only the covering fraction of the denser absorber 
changed significantly and so the model was revised. With the cov¬ 
ering fraction of the denser absorber alone providing the variability, 
the fit statistic stayed approximately the same for an additional free 
parameter and the residuals diminished. The model described the 
data equally well if the column density of the denser absorber was 
allowed to vary with the covering fraction (Fig. [ID. This scenario 
makes sense as two overlapping bodies would alter both the total 
covering fraction and column density along the line of sight. 

A comparative fit was found when the power law parameters 
alone were allowed free to vary between flux states. It was men¬ 
tioned above how this scenario is much more difficult to justify: a 
change in the covering fraction could easily account for short-term 
variability on kilosecond time scales, but it is challenging to imag¬ 
ine a mechanism that would force the slope of the power law to 
fluctuate so quickly. Therefore, the best explanation for the short¬ 
term variability using absorption is a partial covering model where 
a cloud of denser material is moving into the line of sight as another 
cloud with lower density continues to obscure the same solid angle. 

A similar procedure was followed for testing the reflection 
scenario. From past analysis (Section [5.3b it appears the the X-ray 
spectra are power law dominated and that changes in the normal¬ 
ization of the primary emitter itself drive most of the short-term 
variability (Fig.llOb. As the primary emitter brightens or perhaps 
moves along the axis of rotation, the emissivity profile of the disk 
can alter (e.g. Wilkins & Gallo 2015) and material at closer radii 
experiences a change in ionization. When considering the scenario 
of a moving “lamp post” corona emitting isotropically, upward mo¬ 
tion would decrease primary X-ray radiation on the inner disk and 
result in a harder, power law-dominated X-ray spectrum. Down¬ 
ward motion of the corona would cause the disk to be illuminated 
anisotropically due to light bending, increasing the reflected flux 
from the inner region and producing a softer, reflection-dominated 
spectrum. 

The reflection model from Section [53] was applied to the flux- 
resolved spectra (Fig. ED -A good fit could be achieved by allow¬ 
ing the power law component to vary in shape and normalization 
and the blurred reflector to respond with changes in ionization (Ta¬ 
ble |D- Allowing other parameters, such as the reflector normal- 
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Figure 11. The partial covering (top), blurred reflection (middle), and soft 
Comptonization (bottom) models fit to the flux-reduced spectra with their 
respective residuals (data/model) shown below each. The high-flux data are 
black, low-flux are red. The distant reflector, as needed, is shown in blue. 


ization or emissivity profile, to vary did not enhance the fit. Most 
of the change can be attributed to changes in the power law nor¬ 
malization, which seems consistent with the fractional variability 
and difference spectrum analysis. Subtle changes in the ionization 
parameter of the disc arise from changes in the illuminating flux. 


which can explain the low energy excess in the difference spec¬ 
trum. 

Lastly, the soft Comptonization model (Section l5.4b was fit 
to the flux-resolved spectra with only the primary component (F) 
and the radius of the soft Comptonizing region (Rear) allowed to 
vary. Further OPTXAGNF parameters were thawed gradually until 
the most reasonable model was found (Table [7]l. The model that 
gave the lowest xi value allowed Rear, kTe, r, F, and /pi all 
to vary; however, these parameters values were inconsistent with 
those from the average spectral fits. The soft Comptonization model 
that best describes the flux-resolved spectra in a manner consistent 
with the mean spectral fits allowed only Rear, F, and /pi to vary. 
This appears to support the scenario that changes in the primary 
component are responsible for short-term variability in this object: 
as the primary component changes over time, the fraction of the pri¬ 
mary X-rays that are reprocessed in the soft Comptonizing region 
(/pi) also changes, as does the radius of the soft Comptonizing re¬ 
gion (Rcor)- 

The spectrum is dominated by a single component and the am¬ 
plitude and spectral variation are modest in HE 0436-4717, thus 
making more complex timing analysis difficult to justify. That be¬ 
ing said, a lag analysis was performed using three different energy 
bands for comparisons. Given the data quality, lags could only be 
constrained to within a few hundred seconds in each frequency 
band. There was no significant detection of a frequency-dependent 
lag (E. Cackett, private communication), but this could simply be a 
data quality issue. 


7 DISCUSSION 

The multi-epoch, multi-band analysis of HE 0436-4717 has re¬ 
sulted in several scenarios that describe the spectra of the object 
equally well on a purely statistical basis. In the first interpretation, 
emission from the central engine passes through at least one ionized 
body of gas that absorbs the incident radiation. The spectrum pro¬ 
duces an apparent soft-excess that is an artifact of the absorption. 
The second scenario describes the X-ray emission as a composite of 
multiple sources. Here, the primary source from a compact “lamp 
post” corona dominates over emission from blurred, ionized reflec¬ 
tion and a distant, neutral reflector. The Fe Ka line seen in this case 
is a sum of the relativistically broadened iron line originating from 
the inner regions of the accretion disk and a narrow line from the 
distant reflector. The soft-excess is a consequence of blurred ion¬ 
ized reflection from the inner regions of the accretion disk. Lastly, 
it could be that the X-ray spectra of HE 0436-4717 is composed of 
a primary component from a radially extended corona in a “sand¬ 
wich” configuration. The soft-excess in this case is produced by a 
soft Comptonizing region in the inner disk and a distant reflector 
produces the Fe Kq line. 

It is not unusual for these interpretations to model a spectrum 
equally well. Due to the complex nature of AGN, it is reasonable to 
assume that perhaps all three mechanisms (partial covering absorp¬ 
tion, blurred reflection, and soft Comptonization) can be present in 
a single object. However, the degeneracy between the models can 
lead to confusion when attempting to reconstruct source geometry. 
Timing analysis can be key in determining the processes that are 
dominant at a given time. Variability in HE 0436^717 is described 
well with reflection, as it appears that changes in flux of the pri¬ 
mary emitter can be seen influencing both the short- and long-term 
behavior of this object. Assuming an isotropic lamp-post corona, 
changes on kilosecond time scales reflect an increase/decrease in 
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Table 7. Summary of models fit to the flux-resolved spectra. Column 1 lists the individual models and Column 2 the model components. Note that only the 
soft Compton component of the total soft Compton model (Table|^ is shown for simplicity. The model parameters are shown in Column 3 with their values 
for the high- and low-flux states listed in columns 4 and 5, respectively. Parameters linked across the epochs are denoted by dots and those without errors are 
frozen at the given values. 


(1) 

(2) 

(3) 

(4) 

(5) 

Model 

Component 

Parameters 

high-flux state 

low-flux state 

Partial Covering 

Absorber 1 

r 

2.40/“;“ 


Absorption 


3Vh,i (cm-2) 

51.17+^:“? 




CFi 

0.62/“;“ 

U. '^_o.04 


Absorber 2 

3Vh, 2 (cm“2) 

6.00//“ 




CF2 

n Q9+0-06 




Xj/d.o.f. 

1.18/769 


Bluixed Reflection 

primary power law 

r 

2.14±0.01 

2.12 ±0.02 



-^0.1-100 keV 

2.17±0.04 X 10^“ 

1.55 ±0.04 X 10^“ 


blurred reflection 

^in 

5.67 




Ckout 

3.0 




-^br(R'g) 

6.0 




R'ln 

1.79 




sn 

43.5 




Ape (solar) 

0.36 




r 

r primary 




^(ergcm.s“^) 

72 ±4 

53 ±3 



-^0.1-100 keV 

1.12±0.05 X 

0.71 ±0.06 X 10-1’ 


reflection fraction 

7^ 

0.52 ±0.03 

0.45 ± 0.04 


distant Reflector 

Ape (solar) 

1.0 




r 

1.9 




^(ergcm.s“^) 

I.O 




Xj/d.o.f. 

1.12/772 


Soft Comptonization 

soft Compton 

*tBH (M0) 

5.9x10’’ 



component 

Dl(Mpc) 

233.1 




log(L / LEdd) 

-1.05 




logi?out 

3 




a 

0.8 




^?cor ( R'g) 

14.84/j;““ 




kTe (keV) 

0.61+°'^° 




T 

ly >70-1-1.19 

‘■‘<>-1.27 

l-SS/Ha 



r 

2.12±0.01 



/pi 

0.54±0.02 

0.44/“;“ 



xS/d.o.f. 

1.42/772 



primary X-ray flux with respect to the accretion disk. The changes 
in flux could either be intrinsic or due to the motion of the corona 
(e.g. Wilkins et al. 2014; Gallo et al. 2015), which in turn in¬ 
creases/decreases the ionization of the disk at inner radii. Over 
longer, multi-year time scales fluctuations in power law slope can 
be seen which would be caused by physical changes in the corona 
itself (i.e. an alteration in temperature and optical depth). HE 0436- 
4717 remains in a power law-dominated state throughout the ob¬ 
served epochs. When looking at the broader picture of HE 0436- 
4717 and its behavior, the blurred reflection interpretation is the 
most self-consistent model. 

The preference for a reflection scenario in HE 0436^717 
leaves the challenge of explaining the persistently low iron abun¬ 
dance. Observations have shown that most AGN have solar to 
super-solar metallicity values independent of redshift (e.g. Storchi- 
Bergmann et al. 1998, Hamann et al. 2002). In several independent 
studies, it was shown that most AGN are also located in massive 
galaxies with classical bulges (Green & Ho 2004, Barth et al. 2005). 
As stellar mass is strongly correlated with metallicity (Tremonti et 
al. 2004), it is reasonable to expect the majority of AGN to thus 
have solar metallicity values or higher. The question then arises: 
under what circumstances would AGN be observed with (appar¬ 
ently) intrinsic sub-solar metallicities? 

Groves et al. (2006) set out to find such AGN and searched 


SDSS observations of Seyfert 2 galaxies to find those that have sub¬ 
solar metallicities. Out of a sample of around 23,000 they found 
only -^40 candidates based on optical line ratios from the NLR. 
The authors suggest that most AGN are only observed after periods 
of intense star formation in their host galaxies and so the infalling 
gas supplying the central engine is already abundant with heavy 
elements. This implies that AGN with inherently sub-solar metals 
may somehow be feeding off gas that has yet to be processed by 
stars. It could be that HE 0436^717 has an unusually low star 
formation rate for its redshift and/or less Type la supernovae to 
produce iron. Another possibility is that low-metallicity gas from 
outside of the galaxy has made its way to the central black hole 
either through a merger or filament, as in the scenarios proposed 
to explain the unusually low abundances in the NLR of radio-quiet 
quasar HE 2158-0107 (Husemann et al. 2011). 

Other Seyfert 1 galaxies in particular have been known to 
show sub-solar iron abundances, such as NGC 4593 (Reynolds et 
al. 2004), NGC 3227 and 4051 (Patrick et al. 2012), and MCG-h 8- 
11-11 (Bianchi et al. 2010). Patrick et al. (2012) suggest that a per¬ 
sistently low abundance could represent a component of the Fe Ka 
line originating from a Compton thin material such as the BLR or 
NLR. The full-width half-maximum (FWHM) of the iron line was 
calculated for both the pn and merged MOS spectra. The measured 
FWHM corresponded to velocities of ~0.45 c and ~0.31 c for pn 
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and MOS, respectively. GIO measured the FWHM of the optical 
H/3 line in HE 0436-4717 and determined the BLR velocity to be 
around 3990kms~^ or about 0.01 c. Thus, while the exact nature 
of the iron feature in HE 0436-4717 remains uncertain, it is clear 
the line is not originating from the BLR but instead display disk¬ 
like velocities (> 0.1 c, Miller 2007). A narrow iron feature from a 
distant reflector is also required to describe the data in this object. 
The narrow line could originate from the NLR, but it is unlikely 
that such a significant fraction of the radiation originates from the 
region (e.g. Kaspi et al. 2002, Fukazawa et al. 2011). 

Another alternative is the possibility of spallation reducing the 
reflected Fe Ka flux due to bombardment from high-energy (> 
lOMeV) protons. Skibo (1997) explains how abundances of sub- 
Fe metals between 4.5 - 6.4 keV can be enhanced in the nuclear 
regions of AGN at the expense of the Fe line. The author calculates 
the abundances of metals that could fall in the Fe Ka red-wing and 
finds that the corresponding fluxes (specifically those of Ti, V, Cr, 
and Mn) relative to Fe would be large enough to raise their lines 
above the continuum. The combined strength of these sub-Fe lines 
would be comparable to the strength of the Fe line from which these 
lines are produced. If this were the case in HE 0436^717, addi¬ 
tional features in the red-wing should be detected. In particular, Cr 
(at 5.4keV) and Mn (at 5.9 keV) are predicted to have a combined 
flux ~ 40 - 50% that of the observed Fe Ka line - no such features 
are detected in the XMM-Newton spectra. 

Finally, perhaps the intrinsic amount of iron in the central en¬ 
gine is not sub-solar, but the amount of other metals higher since 
the abundance of other metals in the spectral models is fixed at their 
respective solar values. The iron abundance in this object is low 
enough, however, that the other metals in the system would have to 
be super-solar in order to account for the apparent lack of iron. It 
is a challenge to determine a mechanism that would substantially 
increase other metals in a gas without affecting the iron, although 
a lack of Type la supernovae (mentioned previously) may do this. 
More specifically, a decrease in the ratio between Type la/Type II 
supernovae could decrease the relative abundance of iron as both 
produce the heavy metals that decay into iron, however Type la 
more efficiently disperse material into the surrounding interstellar 
medium (i.e. no supernova remnant remains). 

The possibility of HE 0436^717 having intrinsic sub-solar 
iron abundance is intriguing as such objects are rare. Deep analysis 
of the host galaxy would help to further constrain the metallicity of 
this object and also advance understanding of AGN-host relations. 


8 CONCLUSIONS 

A broadband analysis of the Seyfert 1 galaxy HE 0436^717, cov¬ 
ering from optical up to 100 keV is presented. SED modelling sug¬ 
gests the source is unobscured and perhaps its black hole mass 
and distance could be better constrained. A reverberation mapping 
campaign would be advantageous to accurately determine the black 
hole mass in HE 0436^717. 

The multi-epoch X-ray spectra are fitted with Comptonisation, 
blurred reflection, and partial covering models, but do not reveal a 
significantly preferred model. Consideration of the rapid variability 
favors a blurred reflection scenario where the variability is domi¬ 
nated by changes in the power law normalization accompanied by 
corresponding changes in the disk ionisation. However, the blurred 
reflection model indicates a very low iron abundance for HE 0436- 
4717. The low iron abundance does not seem to be a artefact of the 


modelling process, and abundances are not tested in the Compton¬ 
isation or partial covering scenarios. 

Observations of the host galaxy are limited, but would be ben¬ 
eficial to examine any abnormalities in the distribution of abun¬ 
dances in HE 0436^717. Euture X-ray observations with the 
calorimeter Astro-H (Takahashi et al. 2012) could also reveal nar¬ 
row features in the X-ray spectra that would allow abundances to 
be better constrained. In addition to revealing narrow spectral fea¬ 
tures, Astro-H will also allow for simultaneous, broadband soft and 
hard X-ray observations. This should remove some of the ambigu¬ 
ity that is currently present from utilizing BAT data alongside that 
from XMM-Newton . 
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